Coreless and singular vortex lattices in rotating spinor Bose-Einstein condensates 
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We theoretically investigate vortex-lattice phases of rotating spinor Bose-Einstein condensates 
(BEC) with the ferromagnetic spin-interaction by numerically solving the Gross-Pitaevskii equation. 
The spinor BEC under slow rotation can sustain a rich variety of exotic vortices due to the multi- 
component order parameters, such as the Mermin-Ho and Anderson- Toulouse coreless vortices (the 
2-dimensional skyrmion and meron) and the non-axisymmetric vortices with the sifting vortex cores. 
Here, we present the spin texture of various vortex-lattice states at higher rotation rates and in the 
presence of the external magnetic held. In addition, the vortex phase diagram is constructed in the 
plane by the total magnetization M and the external rotation frequency Q by comparing the free 
energies of possible vortices. It is shown that the vortex phase diagram in a M-f2 plane may be 
divided into two categories; (i) the coreless vortex lattice formed by the several types of Mermin-Ho 
vortices and (ii) the vortex lattice filling in the cores with the pure polar (antiferromagnetic) state. 
In particular, it is found that the type-(ii) state forms the composite lattices of coreless and polar- 
core vortices. The difference between the type-(i) and type-(ii) results from the existence of the 
singularity of the spin textures, which may be experimentally confirmed by the spin imaging within 
polarized light recently proposed by Carusotto and Mueller. We also discussed on the stability of 
triangular and square lattice states for rapidly rotating condensates. 

PACS numbers: 03.75.Lm, 03.75.Mn, 67.57.Fg 
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I. INTRODUCTION 



The quantized vortex is one of the hallmark of super- 
fluidity and such the macroscopic quantum phenomenon 
has been studied in many different physical fields ranging 
from the condensed matter physics to neutron stars and 
cosmology. Recently, quantized vortices have been suc- 
cessfully created in Bose-Einstein condensates (BEC's) 
of alkali-atom gases confining in a magnetic trap by sev- 
eral experimental methods 0, 0, IE 13 ■ The static and 
dynamic properties of vortices and vortex-lattices have 
been investigated both theoretically and experimentally 
. Several groups are now able to prepare a vortex array 
with vortices more than one hundred in a BEC 



Further achievement of Bose-Einstein condensation 
was recently done by all optical methods without recourse 
to the magnetic trapping in [10], ^^Rb [Hill El, 
and ^^^Yb atoms 13]. Apart from ^^'^Yb atoms, these 
systems, so-called spinor BEC's, can keep internal de- 
grees of freedom with the hyperfine spin F = 1 or 2. 
The two-body interaction of spin-1 bosons are written as 
dn + gsF ■ F with two parameters gn and gs which de- 
scribe density and spin interactions, respectively |l4lll5j . 
The spin interaction of ^^Na spinor condensates is anti- 
ferromagnetic gs > 8], while it is shown by Klausen 
et al. \&\ that ®^Rb condensates have a ferromagnetic 
nature < 0. Thus, we now have spinor BEC's with 
both ferromagnetic and antiferromagnetic interactions. 

Due to this internal degree of freedom, a rich vari- 
ety of exotic topological excitations have been proposed 



by a large number of authors after the pioneering work 
by Ohmi and Machida Q and Ho 3- monopoles 0, 
[Ta . , 3-dimensional skvrmion s |2C and 2-dimensional 
skyrmions (coreless vortices) ^ 22 , 2^ |3 , alice strings 
(half-quantum vortices) |2£ 



and the other un- 
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conventional vortices |28l |2g|, |30|, Hi Isl Isl Isl l35|. 

Spinor BEC's provide us an opportunity to minutely 
study these exotic properties as a new example of multi- 
component superfluids. Among such exotic state, in par- 
ticular, skyrmion excitations play an important role in 
the other physics systems: the quantum field theory |36j |. 
the quantum Hall system [s^l, superfluid ^He j^sj, ne- 
matic liquid crystals IH, and unconventional supercon- 
ductivity ^3 • The periodic structures of such the topo- 
logical state have been studied in the rotating superfluid 
^He ^38, 41, 44111143, called the Mermin-Ho |4ai and 
Anderson- Toulouse vortices 46] . In the quantum Hall 
system, the textures and stability of skyrmion lattices 
have been also studied m . 

The earlier theoretical study focused on skyrmion exci- 
tations in a spinor BEC has been performed by Khawaja 
and Stoof JO], who pointed out that the 3-dimensional 
skyrmion is not a thermodynamically stable object. The 
numerical analysis of the Gross-Pitaevskii and Bogoli- 
ubov equations under slow rotation, however, led to the 
conclusion that such the topological excitation in the 2- 
dimensional disk is stable and robust in the ferromag- 
netic spin-interaction In addition, at the high ro- 
tation limit, Reijnders et al. ^33] have presented vortex- 
and skyrmion-lattice states for the spin-1 bosons in the 
lowest Landau level. They have also analyzed the ex- 
act ground state near the critical rotation and have pro- 
posed a rich variety of the quantum Hall liquid state. In 
the mean-fleld regime for rapidly rotating spinor BEC's, 
the numerical study of the Ginzburg-Landau equation 
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has been performed within the higher Landau levels [33 ■ 
Here, it has been found that the several types of the vor- 
tex lattice with the shift-core compete with each other, 
which sensitively depends on the spin interaction. 

Previously in our series of papers [2^ |^, we have 
identified that in spinor EEC's with the ferromagnetic 
spin-interaction, the Mermin-Ho vortex is favored under 
slow rotation. This coreless vortex naturally connects to 
the non- vortex state as the total magnetization becomes 
higher, i.e., the scalar condensate limit. Upon increas- 
ing the rotation rate, the stable region of the Mermin-Ho 
vortex shifts toward the low magnetization region, while 
the another vortex state, called the polar-core vortex, 
appears in the high magnetization region. This is the 
vortex state filling the core with the pure polar (anti- 
ferromagnetic) state. In this paper, we will investigate 
the structures of many vortices by numerically solving 
the Gross-Pitaevskii equation. The main purpose of the 
present paper is to address what is the sequence of states 
in a spinor BEC with the higher rotation frequency and 
complete the vortex phase diagram in the plane; the ex- 
ternal rotation frequency f2 versus the magnetization M 
in order to help establishing the properties of textures in 
spinor EEC's. In addition, we will discuss on the stabil- 
ity of two lattice states for rapidly rotating condensates, 
such as the Abrikosov lattice filling the cores with the 
polar state and the square lattice having the continuous 
texture |4ll |. 

This paper is organized as follows. In Sec. II, we first 
present the extended Gross-Pitaevskii equation for spinor 
EEC's, and then explain the numerical procedure to find 
local minima of the energy functional. The spin textures 
and other properties of the favored single- vortex state is 
shown in Sec. HI. In Sec. IV, we present the vortex 
phase diagram in the 51- A/ plane, obtained by comparing 
free energies. We also displayed the detailed structures of 
each ground state. Furthermore, in Sec. V, we show the 
instability of triangular or square lattice states for rapidly 
rotating condensates upon changing the total magnetiza- 
tion (or the external magnetic field) . The conclusion and 
discussions are given in Sec. VI. 



II. THEORETICAL FORMULATION 

A. Gross-Pitaevskii equation 

The Hamiltonian for the F = 1 spinor EEC in a frame 
rotating with the frequency Vl — flz is 0, 0| 



n 



dr 



3 

\ Yl *i {gJjAi + 9sFik ■ Fji} 



ijkl 



(1) 



where "^j and VfJ are the field creation and annihilation 
operators for a boson in the eigenstates of k, I = 

0, ±1). The one-body Hamiltonian is written as 



Mr) 



2to 



v{r) - ni. 



(2) 



with the two-dimensional confinement potential V{r) — 
^moj'^^x^ -\- y"^) and the projection of the angular mo- 
mentum to the z-axis — —ifiixdy — ydx)- The in- 
teraction between atoms with the mass m is character- 
ized by the interaction strengths through the "density" 
channel, c/„ = • °°+^°" , and the "spin" channel, 

= ■ , where oq and 02 are the s-wave scatter- 
ing length with the total spin and 2 channels, respec- 
tively, /i and Bz correspond to the chemical potential 
and magnetic field along the z-axis, respectively. 

The spin angular momentum operators F^ {a = x, y, z) 
with F = l can be expressed in matrices as 



F^ = 



F„ = — I 1 -1 I , (3) 




Fz = 



where the basis is taken as the eigenvector of the spin 
projection along z-axis. These operators satisfy the com- 
mutation relation, [Fa,Ff}] = ieap-yFj. In this basis, 
the field operators for spin-1 bosons are described as 

* = (1'+i,*o,*-l). 

Replacing * in Eq. iQJ by the condensate wavefunction 
■0 = (^) and following the standard procedure, the time- 
dependent Gross-Pitaevskii (CP) equation is obtained as 



d 

i^^V'j = [{h - /-i - jBz 



gnp}Sjk+gs{F) •Fjfe]V'fc,(4) 



where {A) — X^ij V'j*(r)^ijV-'j(r) and the local density 
p{r)=il}^r) ■ ■0(r). 

Here, equilibrium states are found numerically via 
imaginary time propagation of Eq. Q); t —^ t ~ —it, 
starting from arbitrary initial state with random phases 
of each component and with random vortex configura- 
tions. This numerical procedure is equivalent to find the 
local minima of the free energy functional 



1 



dr I (h) 



(5) 



where fi and B^ are interpreted as the Lagrange multi- 
pliers. We use the total number iV2D and the total mag- 
netization M = J2j J drjltpjl"^ as independent variables. 
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Since we assume uniformity along the z direction, the 
order parameter must fulfill the normalization condition 



2D — 



= Y, \ dr|^,(r) 



(6) 



In order to satisfy this condition, the chemical potential /i 
is varied during the numerical iteration. The propagation 
in the imaginary time continues until the fluctuation in 
/i becomes smaller than 10"^" and also that in the norm 
becomes < 10~^ 

The actual calculations are carried out by discretizing 
the two-dimensional space into 100^-400^ mesh. We have 
performed extensive search to find stable vortices, start- 
ing with arbitrary initial vortex configurations for the 
ferromagnetic interaction strength, gs/gn ~ —0.02. In 
addition, we use the following parameters: the mass of 
a ^^Rb atom to=1.44x 10~^^kg, the trapping frequency 
a;/27r=200Hz, and the particle number per unit length 
along the z axis N2t> = 10^/ //m. 



B. Local spin and nematic orders 

It is convenient to introduce a new set of basis, V'a 
{a = x,y^ z) where the quantization axis is taken along the 
a direction [14!|. The transformation from the Cartesian 
representation ipa to ij^j is obtained as 



V'+i \ 


( ^/5 


i 






V^o = 1 


u 





D 




4>-l 


\ ^/2 


i 

%/2 




\4>z 



(7) 



Then it is useful to adopt the following representation 
with real vectors m and n: 



xl) = \fp{m + in). 
We also define the s'pin texture as 
Z = m X n. 



(8) 



(9) 



which describes the direction of the local spin, I = 
{F) / p. Recently, Mueller 'ss'l shows that a tensor Qap = 
'0* (r)-0^(r) can be decomposed as follows: 

Qai3{r)^ieaf:ijl^{r)p{r) + A^Q/3(r), (10) 
N^f,{r) = S^pp{r) - i [{F^F^) + {F^F^)] , (11) 

where Na/3(j) is a symmetric tensor with second-rank 
and describes the spin fluctuation. Furthermore, Caru- 
sotto and Mueller __49] propose that the local values of 
these spin and nematic orders can be imaged by using 
the polarized light. Thus it may be convenient to deflne 
the local nematicity as 



The amplitude of this nematicity is obtained as 
A^(r) = l-|E''W' 



(13) 



which reflects the competing characteristics between the 
local spin order and the local nematic order. 



III. SINGLE- VORTEX STATES 

In the axisymmetric system, the order parameter 
is obtained within the winding number of the j-th 
component Wj and the relative phase aj as ipj — 
\ipj \ exp [i{wj9 + aj)]. From the minimizing of the spin- 
interaction energy, gs{F)'^, it follows that the wind- 
ing and the relative phase satisfy the relation 2wo — 
w+i + W-i and 2ao = a+i -|- a-i -I- mr where n 
is an integer and the even (odd) number corresponds to 
the ferromagnetic (antiferromagnetic) interaction. The 
spinor order parameter then may be written as 



Ipj = y/pjcxp [i{w9 — i{w'6 + «))], 



(14) 



and the spin texture in the x-y plane is given as {lx^ ly) oc 
(cos (w'd + a), sin {w'9 + a)), which is classified with the 
winding number w' and the relative phase a. Here, w 
and w' are related to Wj by Wj = w — jw' , and a and aj 
also satisfies the relation aj — —ja. 



A. Coreless vortices 

In the ferromagnetic gas with ~ 1, m and n in 
Eq. (O are unit vectors with m ± n and hence, these 
three real vectors {l,m,n) form a triad From the 

superfluid velocity v — J2a ["^a^"" ~ naVrUa], the 
local vorticity is obtained as 



V X v = — ^ eafs-yLiyifj) X (V;^), 



(15) 



which implies the Mermin-Ho relation [43. In other 
words, since the original Hamiltonian for the non- 
rotating system with the ferromagnetic spin-interaction 
has the 5*0(3) symmetry the local spins may sweep 
the whole or half unit sphere. The Mermin-Ho (MH) 
vortex is thermodynamically favored under slow rotation 
and the weak magnetic field where the spin interaction 
rather than the contribution of the magnetic field domi- 
nates the free energy, Eq. (0). This coreless vortex takes 
1/7 = 1 and w' = 1, corresponding to the winding com- 
bination (w+i. Wo, W-i) — (0, 1, 2). This is parametrized 
as 



7V= Tr[iV2]/p. 



(12) 



e-'i^'<^+'^) cos2 f 
sin ^ cos ^ 
■ sin^ f 



i(w'0+a) 



(16) 
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where the bending angle f3{r) runs over < /3(r) < tt 
and 9 signifies the polar angle in polar coordinates. The 
spin texture is given as 



sin f3{r) cos {w'O + a) 
— I sin /3(r) sin (^'6* + a) 
cos/3(r) 



(17) 



where f3 has the flexibility for M 23]. As shown in 
Figs.^a), (b), and (c), the MH vortex has a rich variety 
of the two-dimensional type of spin textures by changing 
the relative phases in between spin components: Typical 
textures form the radial disgyration I \\ r (a = 0,7r) and 
the circular disgyration I _L r (a ~ tt/2). The another 
type of the MH vortex is classified as w' ^ —1, i.e., the 
winding combination (2,1,0). As seen in Fig. ^fd), the 
projection of the texture to the x-y plane forms the cross 
disgyration. This texture is also called mixed-twist (MT) 
texture ISIil. These two types (0, 1,2) and (2, 1,0) of 
the coreless vortex completely degenerate at M — 0. 



B. Singular vortices 

In the high magnetization region, the external mag- 
netic field gives the free energy dominate contribution 
beyond the spin interaction energy. The texture then 
depends sensitively on the applied magnetic field. The 
stable vortex in slowly rotating EEC's is the polar-core 
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FIG. 1: The spin textures of isolated Mermin-Ho vortices; (a) 
and (b) radial disgyrations {w' = + 1, a = 0, tt). (c) the circular 
disgyration (ui' = -f 1, Qf = 7r/2). (d) the cross disgyration with 
w' = -1. 
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FIG. 2: The spin amplitude of the /-vector (solid line) and 
local nematicity Miy) (dashed line) of the polar-core vortex. 
The spin texture projected to the x-y plane is shown in the 
inset. 



vortex, which takes the winding it; = and w' = —1, 
corresponding to the combination (1, 0, —I). This vortex 
is thermodynamically favored over two type coreless vor- 
tices in the high magnetization region, which is indepen- 
dent on the strength and the sign of the spin-interaction, 
i.e., the polar-core vortex is the stable object for both 
the ferromagnetic and antiferromagnetic interaction |^ . 
In this configuration, ipQ component with zero winding 
number i^o = occupies the central region of the vortex 
core which is made up by the ipi component. Figure El 
shows the spin amplitude and the local nematicity. This 
state is classified as the winding number w' — —1 and 
thus 2-dimensional texture forms the cross disgyration, 
shown in the inset of Fig.|21 However, this texture differs 
from that of the MT vortex in the foUowings; (i) The 
polar-core vortex has the singularity at the center of the 
cross disgyration where the local nematic order grows up. 
(ii) Since the spin = component is localized in the 
core region whose length scale is the order of the density 
variation characterized by ^„ = h/y/2mnQgn ^ O.l/im, 
the spatial variation of local spin and nematic orders is 
also the order of Here, no denotes the peak density of 
the condensate. As the total magnetization M increases, 
this state continuously connects to the conventional vor- 
tex state in a scalar BEC. 

Recently, Leanhardt et al. 'sSj have created the core- 
less vortex by using the Berry phase engineering method 
[Slf . Bulgakov and Sadreev have also found that the 
polar-core vortex, (1,0,-1), becomes stable in a non- 
rotating harmonic trap with an applied loffe-Pritchard 
magnetic field |52|. 



IV. VORTEX PHASE DIAGRAM 

The vortex phase diagram is calculated by comparing 
the free energy in the plane of the total magnetization 
AI and the external rotation frequency fl. In Fig.|3 we 
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display the resulting phase diagram of the vortex state in 
the M-il plane up to O = 0.4w. The phase diagram for 
slow rotation up to = 0.2a; qualitatively agrees with 
the earlier result in Ref. where a few parameter is 
different from the current system. The rest region of the 
phase diagram can be divided into two categories: (i) 
coreless vortex lattices formed by the MH and/or MT 
textures, and (ii) singular vortex lattices filling the core 
with the polar state. The phase boundary is depicted by 
the solid line in Fig. 13 The phase of continuous vortices 
can be classified as follows: (MH-1) the single MH vortex 
state as seen in Figs.^a)-(c), (MH-2) a pair state of the 
off-centered MH vortex with the radial or circular disgy- 
ration, and (MT-n) the vortex states of n MT vortices 
mixing with some MH vortices (n = 1,2,3). Here we have 
identified the stable phase within the number of the MT 
vortex rather than MH vortex, because in the nonzero M 
region the degeneracy of MH and MT vortices is lifted by 
the presence of the external magnetic field and the role 
of the MT vortex becomes important. The details will 
be discussed in Sec. IV A. The other phases labeled by 
(P-n) in Fig.Ois identified as the singular vortex lattices 
formed by a number of the polar-core vortices. It is noted 
that at the limit of a scalar condensate, M/N = 1, the 
region fl < Q.2uj corresponds to the vortex-free state and 
at the rotation rate fl = QAuj the 4 vortices is energeti- 
cally stable. The composite state of MT and polar-core 
vortices also appear in the narrow M region, labeled by 
(C-1) and (C-2). 

For slowly rotating BEC, fl ^ 0.2w, two axisymmetric 
vortices shown in Figs. and [3 are energetically favored. 
It is seen that the single MH vortex is stable in the low 
magnetization region of Fig. Awhile the single polar-core 
vortex appears in the high magnetization region. The 
separation of this stable region can be qualitatively un- 
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FIG. 3: The vortex phase diagram in the plane by the to- 
tal magnetization M and the external rotation frequency Q,. 
The solid lines represent phase boundaries between different 
vortex states: the non-vortex state, singular vortices, and con- 
tinuous vortices. The details of each phase are noted in the 
text. 



derstood from the energy of the external rotation drive 
—QLz- It is easy to calculate the total angular momen- 
tum Lz of axisymmetric vortices; by using the total num- 
ber N2T) and the total magnetization M, it is simply 
written as 



:W ~ W 



hN2D 



N2B 



(18) 



At M/N2U — 1, MH vortices with w — w' — 1 are equiv- 
alent to the vortex- free state with Lz/fiN2D = 0. As M 
decreases, the angular momentum increases as a linear 
function of M and it reaches to the value of the single 
vortex Lz/fiN2-D = 1 at M/A^2D = 1- Conversely, the 
polar-core vortex with w = and w' = — 1 forms the con- 
ventional vortex with Lz/hN2D = 1 at M/N2U = 1 and 
as M decreases, the angular momentum reaches to the 
zero because of the growth of the spin components with 
the zero and negative winding numbers. This contrast 
behavior for the total magnetization may be reflected by 
the competition of two axisymmetric vortices in the phase 
diagram. For the higher rotation frequencies beyond the 
single vortex region, it is also seen that the stable region 
of singular and coreless vortices are similarly separated 
along the magnetization. The detailed structures of these 
two phases are discussed separately below. 



A. Coreless vortices 

In the low magnetization region, corresponding to the 
zero or weak magnetic field, the energy of the spin inter- 
action is dominated over the magnetization term. Here, 
the important length scale is characterized b y th e spin- 
interaction strength, = h/ y/2mnogs ^ 10^„ which 
results from the competition between the spin interac- 
tion and the kinetic term. This length scale becomes 
the order of spacing between neighboring vortices. This 
requires the continuous spin texture as the energetically 
favored state. Such continuous vortex states can be in- 
terpreted to consist of several MH and MT vortices. 

Figure Ela) shows one example of the stable vortices 
at = 0.35w and M/N2D = where the vector plots are 
the projection of the local spin I to the x-y plane and the 
density map of lz represents the projection to the z-axis. 
This state has the continuous spin texture formed by two 
MH and two MT vortices which are arranged regularly to 
form a square lattice. This spin texture is similar to that 
proposed by Fujita et al. |4lj in connection to superfiuid 
■^He-A phase under rotation. The total density profile 
forms a quite smooth bell shape. 

As M increases, the angular momentum of the MH 
vortex continuously decreases as a function of M while 
that of the MT vortex increases. Therefore, at the finite 
A/, the MT vortex becomes the important object rather 
than the MH vortex. Figures 01a)-(c) show the equilib- 
rium state for the different M: M/N2T, = 0.0 (a), 0.18 
(b), 0.3 (c). At M = 0, four degenerate coreless vor- 
tices forms a square lattice. As M increases, two cores 
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FIG. 4: The textures of continuous vortex lattices near the low magnetization at f2 = 0.35tj: the vortex states with two MH 
and two MT textures at M/N = 0.0 and L^/hN = 2.44 (a), M/N = 0.18 and L^/TiN = 2.47 (b), and at M/N = 0.3 and 
L^/TiN = 2.58 (c). The density map of 1^ is displayed in the left column of (a) - (c) in the range of ~1 (black) and +1 (white). 
The vector plot corresponds to the projection of the Z-vector to the x-y plane, (d)-(f) are the other possible stationary state at 
M/N2U = 0; n = 0.4tj (d), and 0.45aj (e) and (f). All the vortex states have the continuous spin textures and the total density 
profiles always have the smooth bell shape. 



with the MT texture approach to each other. As seen in 
Fig. ^b), the equilibrium state at M/N2D = 0.18 forms 
the double-core of two MT vortices. Here, two MH tex- 
tures with the circular disgyration start to be locked by 
the external magnetic field, which implies that the MH 
vortex continuously deforms into the vortex-free state. 
In Fig. ^c), the equilibrium texture at M/N2T) — 0.3 
and = 0.35a; is displayed has the axisymmetry with 
the higher winding number w = +2 and w' = —2, i.e., 
(4, 2, 0). This can be also regarded as a pair state of the 
MT vortex with w' — —1. The spin texture then can be 
obtained as Z = zcos/3 -I- sin/3(a; cos 26* — ysin26'). Here 
the bending angle j3{r) varies form /3(i?) = around the 



condensate surface to /3(0) = tt at the core. These three 
states are classified as two MT vortices (MT-2) in the 
phase diagram, Fig|2| 



As the rotation rate f2 increases, however, the vortex 
states with various discrete rotational symmetry around 
the BEC center is favored, such as the threefold (b), four- 
fold (c), and fivefold symmetric states (d) in Fig. 21 It 
should be emphasized that all these states have no sin- 
gularities in the spin texture. 
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B. Singular vortices 

As shown in Fig.O the singular vortices labeled by P-n 
are energetically favored over coreless vortices in the high 
magnetization region M/N2t> > 0.5. At the high magneti- 
zation limit M/N2D = Ij all spins are polarized along the 
z-axis. This spin-polarized gas is then described within 
the scalar order parameter and the rotating ground state 
forms the conventional vortex-lattice with the m-fold dis- 
crete rotational symmetry (to = 2,3,4, 5, and 6) l53|. As 
M decreases, the vortex core formed by the spin Fz = +1 
component are filled by the other condensate with the 
spin Fz—0. We depict the density profiles of each con- 
densate with spin Fz = +1 and components a,t fl — OAuj 
and M/N2D = 0.92 in Figs.l^a) and (b). This state con- 
sists of five polar-core vortices, labeled by P-5 in Fig. 13 
The spin = component is then localized in the nar- 
row region which is the order of the length scale of the 
density variation Therefore the sharp peak of the lo- 
cal nematicity appears in the vortex core region, shown 
in Fig. |Sfc). As seen in Fig. |SJd), the almost spins are 
locked by the external magnetic field and the spin tex- 
ture can continuously vary around the cores. This forms 
the cross disgyration having the singularities at the cen- 
ter of the disgyration. The length scale is characterized 
by ^ and thus the spatial variation of the spin tex- 
ture is much shorter than the spacing between vortices. 
This results in that there is no correlation between each 
polar-core vortex. 

Figures lUe) and (f) show the density profiles of each 
component, p+i and po, at Q — OAoj and M/N2D — 0.71. 
As M decreases, the spin component po spreads in be- 
tween vortices formed by the spin Fz = +1 component. 
This leads to the growth of the local nematicity, shown 
in Fig. [S{g). The spin texture is depicted in Fig. Eth) 
where the spin texture forms five cross disgyrations hav- 
ing the broader spatial variation than that in Fig.[3Jd). 
Then, the coreless vortex with the radial disgyration is 
spontaneously created in the center of five polar-core vor- 
tices in order to smoothly connect the five cross disgy- 
rations. This is similar to the texture of the coreless 
vortices shown in Fig. 0] (f). However the length scale 
of the spatial variation of the texture is characterized by 
the order of the density variation ^„ rather than ^s- In 
addition, the point singularities exist in the spin texture. 
It should be noted that the system keeps the discrete 
rotational symmetry even if the total magnetization M 
further decreases, which is the common feature to the 
other ground states labeled by P-n in Fig. |3] (rt = 2, 3, 
and 4). 



V. VORTEX LATTICES 

Here we show one example of the equilibrium state for 
rapidly rotating EEC's. Under slow rotation, the ground 
state at M/N2T1 ~ consists of coreless MH and/or MT 
vortices as seen in Sec. IV. The initial state is taken 



as the continuous texture periodically arranged to the 
square lattice, whose unit cell consists of two MH and 
two MT vortices as seen in Fig. EJa). By numerically 
solving the GP equation at = CQw in the absence of 
the external magnetic field, the square lattice is given as 
the equilibrium state, shown in Fig.EJa). This lattice is 
constructed from two sublattices of the MH vortex with 
the circular disgyration and the MT vortex with the cross 
disgyration, and the local spins on the two sublattice 
sites are locked to the alternative direction z and — i, 
respectively. It is however found in the equilibrium spin 
that there is slight distortion from a square array around 
the condensate surface where the Thomas-Fermi radius 
is ~ 5.8/xm. 

In the presence of the weak magnetic field, the same 
initial state yields the another equilibrium state. The 
equilibrium texture and the total density profile at 
M/N2D = 0.3 is depicted in Fig. E^b) and (c), respec- 
tively. At the EEC center, there is the non-singular 
skyrmion texture with the higher winding number whose 
spin texture is characterized by Eq. (|17|l with w = +2 
and w' = +2. This is also identified as a pair of the MH 
vortex with the circular disgyration. In addition, two 
MT vortices with the cross disgyration form the double- 
core state and are located around the center. It is how- 
ever found that the other vortices with the cross disgy- 
ration appear around the condensate surface. They have 
point singularities and thus the local nematicity shown in 
Fig. Eld) grows up at the cores. This vortex lattice can 
be identified as the composite lattice of the coreless vor- 
tex and the polar-core vortex. It is noted that two length 
scales of the spin variation exist in this lattice sate, ^„ 
and ^s- 

On the other hand, we discuss on the stability of the 
vortex lattice in the high M region for rapidly rotating 
EEC's. Figures [7|[a)-(d) show the structure of the equi- 
librium state at n — 0.9lu and M/N2D — 0.95: the total 
density profiles (a), the spin texture (b) and (c), and the 
local nematicity (d). The initial state in the calculation 
of the GP equation is taken as the vortex configuration 
with the hexagonal symmetry which is the most favored 
configuration in rotating scalar EEC's, corresponding to 
the high magnetization limit M/N — 1. The resulting 
equilibrium state consists of the vortices filling the cores 
with the polar state and keeps the hexagonal symmetry. 
As shown in Fig. Eld) , the local nematic order grows up 
in the narrow core region where the spin texture has sin- 
gularities. This length scale is characterized by the order 
of the density variation 

Figures [7Je) shows the total density profiles at Vl ~ 
0.9w and M/N2T1 = 0.59. It is seen that the symme- 
try of the lattice drastically changes from the hexagonal 
symmetry to square one, where the initial state is taken 
as the vortex lattice regularly arranged to the hexagonal 
symmetry. Such the transition occurs at the total mag- 
netization M /N2T> ~ 0.93 at $7 = 0.9ti;. The spin texture 
projected to the z-axis is depicted in Fig. [T^f ) and (g). 
As seen in Fig. [7fg), the coreless MH vortex with the 
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FIG. 5; The equilibrium state at ^ = OAlo and M/N2-D = 0.92: The density profiles of the spin = +1 (a) and = (b) 
components, the local nematicity J\f (c) in the region — 6^m < x,y < +6/im. Figure (d) shows the vector in the central region 
— 4/im < x,y < +4^m. In the left figure, the density map of Iz is shown in the range of (black) and +1 (white), and the 
vector plot in the right one corresponds to the projection of the Z- vector to the x-y plane, (e)-(h) correspond to the equilibrium 
state at = 0.4tj and M/N20 = 0.71. The local nematicity in (c) and (g) is displayed in the range of 0.5 (black) and 1 (white). 



circular disgyration is spontaneously created in spacing 
between polar-core vortices with the cross disgyration. 
In addition, the spatial distribution of the local nematic- 
ity spreads over spacing between neighboring vortices, 
shown in Fig. |2Jh). These result from the spatial conti- 
nuity of the spin texture, i.e., the ferromagnetic feature 
of the spin interaction. Such the drastic change of the 
symmetry of the vortex lattice has been also reported in 
two-component EEC's [ssHSGj . 

VI. CONCLUSIONS 

In this paper, we have presented the results of numer- 
ical calculations of the vortex-lattice states in rotating 
EEC's with the ferromagnetic interaction. In the ab- 
sence or presence of the external magnetic field, we have 
studied the local properties of the stable states, such as 
the spin texture and the local nematicity. In addition, 
the vortex phase diagram has been constructed in the 
plane by the external rotation frequency SI and the total 
magnetization M . The stable phase can be divided into 
two categories: (i) coreless vortices with the non-singular 
spin texture and (ii) singular lattices with polar (antifer- 
romagnetic) cores. In the resulting phase diagram, the 
competition between these two phases has been found. 

Coreless lattices which are the most favored state 
in the zero magnetic field are formed by several types 
of the 2-dimensional disgyration, corresponding to 2- 



dimensional skyrmion lattices. In addition, for rapidly 
rotating EEC's, the coreless vortex regularly arranged to 
the square lattice has been presented and the stability has 
been discussed. This square lattice becomes unstable for 
applying the weak magnetic field and then, neighboring 
MT vortices are paired and form the double-core lattice 
or the skyrmion lattice with the higher winding number. 

In the high M region, it has been also demonstrated 
that the singular lattices are energetically favored over 
coreless vortices. The projection of the texture to x-y 
plane in the singular vortex is the same as that in the 
coreless MT vortex. With decreasing M, the spatial vari- 
ation of the spin texture becomes broader than the order 
of the density variation ^„ and the correlation between 
neighboring vortices plays an essential role. As the re- 
sult, the drastic change of the hexagonal lattice into the 
square one has been demonstrated, which may be ex- 
perimentally identified by using the spin imaging within 
polarized light proposed by Carusotto and Mueller [j^. 
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